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Abstract: MINOS is the first large magnetic detector deep underground and is the first to measure 
the muon charge ratio with high statistics in the region near 1 TeV.[l] An approximate formula for the 
muon charge ratio can be expressed in terms of = 1 15 GeV, ex = 850 GeV and ^^ur/ace q jj^^ 
implications for K production in the atmosphere will be discussed. 
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Introduction 



Defining A 



A = 



0.14R72-7 

M 

cni^ s sr GeV 



(1) 



the muon energy spectrum underground can be ex- 
pressed as [2]: 
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[dN] 

(2) 

where the two terms give the contributions of 
charged pions and kaons; the values of e = 
mcho /t are close to the energies where the prob- 
ability of meson interaction and decay are equal: 
= 115 GeV and = 850 GeV. The zenith an- 
gle is 6 and in experiments with a flat overburden, 
the largest muon intensity comes at cos^? = 1. It 
is possible to generalize this equation to separately 
consider /i+ and /^^ by introducing constants /^r 
and Jk- 
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The measured charge ratio in a bin of surface muon 
energy and zenith angle can then be expressed as 




(3) 



115 GcV 



850 GcV 



An interesting feature of Equation 3 is that it de- 
pends only on the product ij^ur/ace q ^j^^ ^^^^ 
otherwise on E or 0. This combination of terms 
controls the relative portions of interaction and de- 
cay for both tt's and K's. Thus at a fixed value of 
j^surface Q jj^g j-^j-Jq q£ from tt's and K's 

is constant. Equation 3 postulates an energy in- 
dependent 7r+/7r^ ratio related to an energy 
independent / ratio related to fx and an 
energy independent tt/K ratio embodied in the 
0.054. These are all reasonable as a consequence 
of Feynman scaling. (At energies above 10 TeV, 
charm and the changing cosmic ray chemical com- 
position could also affect the muon charge ratio.) 
In this simple but powerful parameterization, the 
charge ratio depends only on the relative number of 
tt's and K's that contribute to muons measured un- 
derground. The relative contribution to the muon 
intensity from tt's and K's is affected by the rel- 
ative probability of hadron interaction and decay, 
which vary with both energy and angle. It varies 
as a function of energy because of the hadron life- 
times, and it varies as a function of angle because 
of the interaction probability in the atmosphere. 
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But as can be seen from inspection of Equation 3, 
at any fixed value of £:s«'"/«ce ^jjg ^^^^ 
muons from tt and K decay is also fixed, even if E 
or 6 is varied. 

Measurement of "'^^"'^^ cos 6* under- 
ground 

The energy of a muon underground can be re- 
lated to the energy at the surface by Egurface = 
Eunderground+Eioss The energy loss for muons in 
MINOS has been calculated using the energy loss 
equation, dE/dx = a{E) + h{E) x E. Many high 
energy muons which do reach the detector do not 
bend enough in MINOS' magnetic field to measure 
the curvature and hence charge. This can be quan- 
tified in terms of the maximum detectable momen- 
tum, which depends on the geometry of the track 
through the detector [3]. A magnetic detector un- 
derground can measure the charge of muons that 
have the lowest energies Eunderground when they 
reach the detector. 

MINOS has made measurements of the charge ra- 
tio in the far detector [1] which is 710 meters un- 
derground or 2070 m.w.e., and also has yielded a 
preliminary result using its near detector [4] which 
is 100 m underground or 225 m.w.e. (Note that the 
far detector is under a small ridge.) 

As a crude approximation, we can consider the 
case when: 1) the parameter a is constant, 2) the 
parameter b is zero. 3) the rock density is con- 
stant in all directions, 4) the surface is flat and 
5) the maximum detectable momentum is negligi- 
ble compared to the energy loss from the surface. 
The first four assumptions are equivalent to saying 
that the energy loss is proportional to the overbur- 
den which only depends on zenith angle. In that 
case Eioss = Emin/ cos 9 where Emin is the min- 
imum energy loss for a vertical cosmic ray muon. 
The fifth assumption implies that the underground 
energy of muons used for the charge ratio mea- 
surement in MINOS was much smaller than their 
surface momentum. The MINOS maximum de- 
tectable momentum, defined at one sigma, is dis- 
cussed in these proceedings [3] and a cut was cho- 
sen at 2.2 sigma. 

In the limit corresponding to these assumptions, 
the frequency distribution in £;«"'"/«ce ^.Qg q ^^^^ 



we would measure is a delta function at Emin- At 
large zenith angles, the surface energy increases 
due to the cos 9 dependence of the slant depth, but 

j^surface Q COUStaUt. 

This can be compared with the measured E and 
^surface ^^^Q distributions from MINOS for the 
far detector in Figure 1 and for a preliminary re- 
sult from the near detector in Figure 2. Both 
^surface ^gg Q distributions are much narrower 
than the corresponding Egurface distributions. The 
largest contribution to the width of the measured 
^surface ^^g Q distribution for the far detector is 
the b(E) term in the energy loss, while the largest 
contribution to the much narrower near detector 
distribution is the larger ratio of maximum de- 
tectable momentum to energy loss in the overbur- 
den. 

Implication for Kaons 

The charge ratio of muons that we measure from 
TT decay is r^r =t"%- and from kaon decay fx 

= ,-^'i ■ The value 0.054 is related to the relative 

kaon/pion mtensity in Gaisser's model. [2] 

A fit was made to MINOS data and the published 
L3-hC ratio data [5] to this parameterization for r. 
A chi-square fit was done in 178 bins of 
vs. cos 9. The fit gave U = 0.555 ± 0.002 and 
fx = 0.667 ± 0.007. The errors in and fx are 
correlated. [1] Using these numbers, the charge 
ratio versus "^^e ^Qg^ plotted in Figure 3. 
These values for /^r and fx imply that the muon 
charge ratio from pion decay is r = 1.25 and that 
the muon charge ratio from kaon decay is about 
r = 2.0. 

As noted above, there is some £;f "''/«ce ^^^g q 
ation in the MINOS far detector data at a value of 
j^surface Q where the curve in Figure 3 is vary- 
ing. We have performed a fit solely to the MINOS 
far detector data, and obtain f^, = 0.549 ± 0.012 
and fx = 0.702 ± 0.049, which imply = 1.22 
and rx = 2.36. The errors on and fx in this 
fit are almost completely correlated and are shown 
in Figure 4. A fit which also incorporates the pre- 
liminary MINOS near detector charge ratio will be 
presented at the 30th ICRC. 

Equation 3 has asymptotic values at both low 
energy and high energy. Note that both val- 
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ues have contributions from both tt and K de- 
cays. For the above MINOS-only fit, the low 
^sttr/ace ^,yg g charge ratio is r = 1.27 and the 
high ij^tir/ace ^ charge ratio is r = 1.45. The 
parameterization does contain enough physics to 
(so far) satisfactorily represent the world's under- 
ground muon charge ratio data. Furthermore, it is a 
reasonable way to parameterize the energy depen- 
dency of the charge ratio, and seems preferable to 
a Unear fit in E or log E. We suggest that the new 
MINOS data, equation 3, and the fitted values for 
the model parameters near 1 TeV will be useful in 
the next round of atmospheric neutrino modeling 
in this energy range. It is also now apparent that 
future precise muon charge ratio data in the 0.2- 
0.6 TeV energy range will be useful for verifying 
the £;™''/<"=e Q dependency of the charge ratio. 

Our fits to r,r give values near expectation [6]. Our 
fit to the rK = K^/K~ charge ratio in atmo- 
spheric showers yield values just above 2. It is 
clearly difficult to directly measure this ratio. Eq. 3 
provides a direction for future study of this subject. 
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Figure 1: Disti-ibution of E and cos 6* for 

MINOS data muons in the far detector, after cuts. 
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Figure 2: Preliminary distribution of E and 
^surface Q MINOS data muons in the near 
detector, after cuts. 
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Figure 3: Behavior of Equation 3 for •^"^'^ cos 6* from 1 GeV to 10 TeV. The critical energies of tt's 
and K's are shown. The regions where the MINOS near and far detectors are sensitive are also indicated. 
Values of /^r and /k from the MINOS/L3+C fit were used. 




Figure 4: Distribution of surface in the /^r - fx plane, showing the error of the MINOS only fit and the 
correlation between and fx ■ 
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